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The Continuum Discretized Coupled Channels (CDCC) method is a well established theory for 
direct nuclear reactions which includes breakup to all orders. Alternatively, the 3-body problem 
can be solved exactly within the Faddeev formalism which explicitly includes breakup and transfer 
channels to all orders. With the aim to understand how CDCC compares with the exact 3-body 
Faddeev formulation, we study deuteron induced reactions on: i) 10 Be at Ed = 21.4,40.9 and 71 
MeV; ii) 12 C at E d = 12 and 56 MeV; and iii) 48 Ca at E d = 56 MeV. We calculate elastic, transfer 
and breakup cross sections. Overall, the discrepancies found for elastic scattering are small with the 
exception of very backward angles. For transfer cross sections at low energy ~10 MeV/u, CDCC is in 
good agreement with the Faddeev-type results and the discrepancy increases with beam energy. On 
the contrary, breakup observables obtained with CDCC are in good agreement with Faddeev-type 
results for all but the lower energies considered here. 

PACS numbers: 24.10.Ht; 24.10.Eq; 25. 55. Hp 
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The construction and upgrade of various radioactive 
beam facilities worldwide continue to open up new av- 
enues for the study of direct nuclear reactions involving 
rare isotopes H 0. At the low energy regime (around 
10-20 MeV/u) one can explore the structure of the ex- 
otic nucleus in detail. In this context, deuteron induced 
reactions in inverse kinematics are of great interest be- 
cause of the reduced Coulomb barrier, when compared to 
heavier targets. These include, elastic, inelastic, trans- 
fer and breakup reactions. Single neutron transfer using 
(p,d) or (d,p) is also becoming a common tool to explore 
neutron capture for astrophysics (e.g. Q). 

Given the long history of using deuteron induced reac- 
tions (in normal kinematics) for studying nuclear struc- 
ture, a number of reaction theories are well established. 
With the understanding that breakup strongly affects the 
reaction mechanism when the loosely bound deuteron is 
involved-theories including deuteron breakup were devel- 
oped [J-jy]. However, until recently the accuracy of the 
methods has not been quantified. Because the standard 
reaction methods are not based on any normal ordering 
perturbation theory, it is not possible to estimate the ac- 
curacy intrinsically. One relies on an independent exact 
formulation of the problem. 

Faddeev devised a method which couples rearrange- 
ment, elastic and breakup channels Q by expanding the 
three-body wavefunction into three components. For any 
given three-body Hamiltonian, the Faddeev method pro- 
vides the exact solution. For few-nucleon reactions an 
integral equation formulation for the transition opera- 
tors @, Q proposed by Alt, Grassberger, and Sandhas 
(AGS) is commonly used in the momentum-space frame- 
work. One critical issue in this method is the well known 



mques |1U|, [ll|j have been developed to handle Coulomb, 
which opened the path for applications to deuteron in- 
duced nuclear reactions [l2l . [l3j . This method can now 
be used to benchmark existing nuclear reaction theories. 

One method of practical use in (d,p) and (p,d) reac- 
tions, is the finite range adiabatic method (ADWA) Q 
which takes into account deuteron breakup to all orders 
in the transfer channel while making an important sim- 
plification of the continuum. ADWA has been studied 
in detail in 14 and compared to exact Faddeev predic- 
tions in 3, 16 . Those studies demonstrate the validity 
of ADWA for (d,p) reactions at low energies (E ~ 10 
MeV/u), populating states with low an gula r momentum 
and small separation energy. Results in [16J indicate that 
the accuracy of ADWA deteriorates with increasing beam 
energy, angular momentum and separation energy of the 
final state. 

ADWA is the precursor of the Continuum Discretized 
Coupled Channel Method (CDCC) @], which docs not 
make approximations in the treatment of the breakup 
states. While ADWA provides a method to compute 
transfer alone, CDCC can be used for all, elastic, breakup 
and transfer. CDCC is now a standard tool to study 
deuteron induced breakup reactions [l7|. In CDCC, the 
3-body scattering problem is solved through a set of cou- 
pled Schrddingcr-likc differential equations. The breakup 
is included to all orders by expanding the full three-body 
wave function in terms of a complete basis of the projec- 
tile's bound and continuum states. A discretization by 
averaging over energy is common to handle the coupling 
between the continuum states. In addition to elastic and 
breakup, attempts have been made to study transfer re- 
actions using the CDCC wavefunction [12, Il8l - [2l| . In 
these methods, the exact three-body wavefunction in the 
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transfer amplitude is replaced by the CDCC wavefunc- 
tion. Thus, starting from one single three-body Hamilto- 
nian, CDCC can provide, in a consistent manner, elastic, 
breakup and transfer cross sections. 

In theory, CDCC should provide an exact solution to 
the three-body problem. In practice, due to the trunca- 
tion of the model space this may not be true. For exam- 
ple, it is easy to demonstrate that the truncated deuteron 
CDCC wavefunction does not provide the correct asymp- 
totic form for a bound state in a rearrangement channel. 
A first attempt to compare CDCC and Faddeev was per- 
formed in Ref. [l2[ . Good agreement was found for elas- 
tic scattering and breakup for deuterons on 12 C and 58 Ni 
while disagreement was observed for 11 Be on protons at 
E = 38 MeV/u in both breakup and transfer observ- 
ables. The reaction calculations for 12 C and 58 Ni did 
not include bound states in the rearrangement channels 
and therefore one may naively expect the CDCC solution 
to provide the correct answer to the problem. The dis- 
agreement found for the single reaction where a bound 
transfer state was included, questioned the reliability of 
CDCC and called for a more systematic study. 

The aim of this work is to quantify the accuracy of 
CDCC in computing elastic, breakup and transfer cross 
sections and establish a range of validity. This will be 
done by direct comparison with exact Faddeev- AGS cal- 
culations (FAGS). Our test cases include deuteron scat- 
tering from 10 Be, 12 C, and 48 Ca. These reactions are 
studied as a function of beam energy, chosen to match 
available experimental data. 

In Section II, we describe the details of the CDCC and 
AGS formalisms. In Section III we present the numer- 
ical details of the calculations, followed by the results 
and discussions, in Section IV. Finally, the summary and 
conclusions arc drawn in Section V. 



A. The CDCC Method 

In the CDCC approach, using the above mentioned 3- 
body Hamiltonian (Eq. [TJ, the 3-body wave equation of 
the total system is written as 
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where r and R are the Jacobi coordinates for the p-n and 
d-A systems, respectively. Here E is the total energy of 
the system in the center of mass (cm.) frame. 

Keeping in mind that the deuteron is a loosely bound 
system, we consider deuteron breakup states explicitly. 
The eigenstates of the deuteron, <fi(r) satisfy the eigen- 
value equation (T r + V np ) (f>(r) = ecf>(r). To reduce the 
computational cost, in this work we ignore the spins of 
the particles. We perform the standard partial wave de- 
composition in terms of the orbital angular momentum of 
neutron relative to the proton (I) and the discretization of 
the deuteron continuum (indexed i) by slicing the contin- 
uum into bins and averaging over momentum [|| [l7l |22| . 
In this manner, the projectile is represented by: 



■M = E 



i Uji(r) , , 



(3) 



where i = refers to the deuteron ground state and i > 
to the np continuum bin i corresponding to momentum 
ki = y/ 2 fi pn t% I H 2 with jjL pn being the reduced mass of 
proton and neutron. The projectile continuum is trun- 
cated by restricting the orbital momentum (I < l ma x) 
and the energy of the np breakup bins (e$ < e max ). The 
limits are chosen such that the observables of interest are 
converged. 

A partial wave decomposition in terms of the relative 
angular momentum between projectile and target (L) is 
also performed. Thus, the full 3-body CDCC wave func- 
tion can be written as: 



II. FORMALISM 

For describing the A(d,d)A, A(d,p)B, and A(d,pn)A 
reactions we consider an effective 3-body problem p + 
n + A , such that B is the ground state B = n + A. The 
Hamiltonian is then written as 



1> 



CDCC 



(r,R) =E 



i l+L ^Xa(R)Y la (r)Y La (K) 



(4) 

where we use a as a general index representing the quan- 
tum numbers a = {ilL}. 

Substituting Eq.Q in Eq.®, a set of coupled-channel 
equations is obtained for Xa(R), 
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where the kinetic energy operators, Tr and T r can be ex- 
pressed in any of the three Jacobi coordinate pairs (22|. 
The deuteron is modeled by a real potential, V pn , repro- 
ducing its binding energy and the np low energy phase 
shifts. The p-A and n-A interactions are phenomenolog- 
ical optical potentials, U p a and t/ n A, respectively. They 
include absorption from all the channels which are oth- 
erwise not included in the model space. 



h 2 f d 2 _ L(L+1) 
2ndA \dR 2 R 2 



+ V aa {R) + e i -E 



+ E %L '~ L V *cc-{R)Xcx'{R) = ,(5) 

where /idA is the reduced mass of d+A, and E a = E — ei. 
Eq.© is solved with the scattering boundary conditions 
for large R: 



Xa (R) -> */2 [H-{KR)5 aaz - H+(KR)S aai ] 



(6) 
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Here, H~ and H + are the Coulomb Hankel function [22[ 
and S aai is the S-matrix. The coupling potentials in 
Eq.([5]) arc defined by: 



V aa ,(R) = (Mr)\U pA + U nA \<f> llv (r)) 



(7) 



and contain both nuclear and Coulomb parts, which are 
expanded in terms of multipoles Q up to Q ma x ■ 

As mentioned before, the CDCC wavefunction Eq.Q 
can be used in the exact transition amplitude for 
A(d,p)B. The post form for the transition amplitude 
is then: 



T 



\ X P B ^ nA I V P n + U P A 



CDCC 



), (8) 



where x p b * s ^ ne distorted wave generated by the aux- 
iliary potential U p b (depending on R', the p-B relative 
coordinate), and </> n A is the bound state for the final n-A 
system. The difference (U p a — U p b) is called the remnant 
term, which is sometimes neglected for convenience. It 
should be stressed that, in the case of exotic nuclei, it can 
be very significant. In our study it is always included. 



B. The Faddeev-AGS method 

The Faddecv formalism explicitly includes the 
deuteron d + A, the proton (nA) + p, and the neutron 
(pA) + n components. Three sets of Jacobi variables are 
used such that each component is treated in its proper 
basis. In the AGS approach the three particle scatter- 
ing is described in terms of the transition operators Tp a , 
where a (j3) corresponds to the initial (final) channel in 
the odd-man-out notation. The transition operators obey 
the AGS equations Q that are a system of coupled inte- 
gral equations 



T, 
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Here 5p a = 1 — 6p a and Go = (E + iO — Hq)^ 1 is the free 
resolvent with E being the total energy in the three-body 
cm. system. The two-body transition operator for each 
interacting pair is derived from the pair potential v 7 via 
the Lippmann-Schwinger equation 
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The scattering amplitude Xp a is the on-shell matrix ele- 
ment of Tp a calculated between initial and final channel 
states as 



X-Pa = {4>p\ T l3a\(t>a 



(11) 



In the case of elastic scattering j3 = a, while (3 = for 
breakup and ^ (3 ^ a for transfer reactions. The AGS 
equations §§§ are solved in momentum space using par- 
tial wave decomposition. The proton-nucleus Coulomb 
interaction is treated using the method of screening and 
renormalization fiol 111. Further numerical details can 



be found in Rcfs. 23 - 26 



III. NUMERICAL DETAILS AND INPUTS 

In this section we present the pair interactions nec- 
essary to define the three-body Hamiltonian as well as 
some details for the model space used. Both CDCC and 
Faddeev-AGS are computationally expensive, thus, with- 
out sacrificing the final goal, we make simplifications to 
the p + n + A system. To reduce the number of channels, 
we have neglected the spins of the particles and thus all 
the spin-orbit terms in the interactions are not included. 
Also, a single gaussian interaction is used to describe the 
deuteron and its continuum 27 



V(r) = -V e- (r/ro)2 



(12) 



with V = 72.15 MeV and r = 1.484 fm. The poten- 
tial reproduces the binding energy of deuteron and low- 
energy p-n 3 Si phase shifts but fails heavily in other par- 
tial waves. While the comparison of CDCC and Faddeev- 
AGS is meaningful given that the same Hamiltonian is 
used in both formulations, a direct comparison of the pre- 
dicted cross sections with data should be avoided. We 
have chosen for all our test cases energies at which data 
exists and thus this study retains its relevance to exper- 
iments. 

An important ingredient to the three-body Hamilto- 
nian, the nucleon-nucleus optical potentials, are all taken 
from the global parameterization CH89 [28| . The optical 
potentials are energy dependent and should be considered 
within the context of the reaction model. In CDCC calcu- 
lations, U p a and U n A are calculated at half the deuteron 
incident energy Ed/2, while the auxiliary potential, C/ p b, 
is calculated at the proton energy E p corresponding to 
the exit channel. The Faddeev-AGS calculations with 
exactly the same choice of U p a and U n A as in CDCC 
yield results for elastic and breakup observables only as 
there is no transfer channel with complex U n A ■ We there- 
fore perform also the Faddeev-AGS calculations where a 
real neutron-nucleus potential is used in the partial wave 
corresponding to the final bound state. The same in- 
teraction is used to produce the neutron-nucleus bound 
state cf> n A in Eq.©. We take a Woods-Saxon central 
potential with standard geometry: radius tq = 1.25 fm 
and diffuseness a = 0.65 fm. The depths are adjusted to 
reproduce the neutron separation energy of the relevant 
state. Details are given in Table |U 

For both CDCC and Faddeev-AGS, convergence needs 
to be studied in detail. While elastic and transfer observ- 
ables offered a minor challenge from the point of view of 
convergence, the same cannot be said of breakup dis- 
tributions. Even though spins were neglected to reduce 
computational cost, breakup calculations proved to be 
at the limit of our numerical capabilities. Due to techni- 
cal difficulties in the Faddeev-AGS, breakup calculations 
were performed without Coulomb. In Section IIII Al we 
present the CDCC model space based on transfer ob- 
servables. For those cases in which the model space was 
insufficient to ensure convergence of the breakup cross 
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Nucleus (A) nl S n (MeV) V nA (MeV) 



J Be 



^Ca 



2.s 



2p 



0.504 



4.947 



5.146 



57.064 



39.547 



48.905 



TABLE I: Parameters of the n- A binding potential used in the 
CDCC calculations. The Coulomb radius and the diffuseness 
are taken to be ra = 1.25 fm and ao = 0.65 fm, respectively. S n 
is the neutron separation energy, V n A is the depth of the n-A 
binding potential and nl are the quantum numbers describing 
the ground state of the nucleus B=A+n. 



sections, the predictions are presented with bands instead 
of lines, estimated on the convergence behavior with in- 
creasing I. In Section UlIBI we present the Faddeev-AGS 
model space chosen to ensure convergence of all calcu- 
lated observables. 



A. CDCC Model space 

The 10 Be(d,p) 11 Be reaction is studied at three beam 
energies E d = 21.4, 40.9 and 71 MeV. The calculations 
include p-n partial waves up to l max = 4 and Qmax = 4 
multipolcs in the expansion of the coupling potentials. 
The coupled equations are integrated up to R ma x = 60 fm 
with total angular momentum J ma x = 40 for E c i = 21.4 
MeV and J max = 50 for E d = 40.9 and 71 MeV. 

Calculations for 12 C(d,p) 13 C reaction includes p-n par- 
tial waves I < 6 and coupling potentials expanded up to 
Qmax = 6. The coupled equations arc integrated up to 
Rmax = 60 fm with J max = 30 for E c i = 12 MeV and 
J max = 50 for E d = 56 MeV. 

Calculations for 48 Ca(d,p) 49 Ca reaction includes p-n 
partial waves I < 6 and coupling potentials expanded up 
to Qmax = 6. The coupled equations are integrated up 
to R max = 80 fm with J max = 70 for E d = 56 MeV. 

For all calculations, the radial integration was per- 
formed in steps of 0.06 fm and the momentum average 
defining the bins was performed on a mesh with 20 points. 
The schematic detail of the bin structure used for the 
discretized deuteron continuum in these calculations, is 
shown in Fig. [T] 

When considering transfer, in addition to the above 
mentioned model space, the non-locality of the trans- 
fer kernel needs to be considered. The details of the 
non-local parameters used are presented in Table [TT1 Due 
to the strong repulsive couplings at short distances, for 
10 Bc(d,p) 11 Be at E^ = 71 MeV, it was necessary to in- 
troduce an L-dcpcndent internal cut in the couplings. 
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FIG. 1: Schematic detail of deuteron continuum used in the 
calculations, showing number of bins and the energy range as 
a function of I: (a) d+ 10 Be, (b) d+ 12 C and (c) d+ 48 Ca. 



B. Faddeev-AGS Model space 

The AGS equations are solved using partial wave de- 
composition. However, in contrast to CDCC, three bases 
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target E d (MeV) Centration (fm) Non-local width (fm) 





21.4 


2.16 


14 


10 Be 


40.9 


1.92 


20 




71 


1.32 


12 


12 C 


12 


-0.96 


6 




56 


-0.66 


7 


48 Ca 


56 


-0.78 


7 



TABLE II: Parameters of the non-local transfer kernel used 
in the CDCC calculations (details can be found in [2!|). 



corresponding to the three choices of Jacobi momenta are 
used. The maximum we include for the p-n partial waves 
is Imax = 5. For the 10 Be target, n- 10 Be partial waves 
up to Imax = 8 and p- 10 Be partial waves up to l max = 16 
are included. For the 12 C target, n- 12 C partial waves up 
to l m ax = 8 and p- 12 C partial waves up to l max — 18 are 
included. For the 48 Ca target, n- 48 Ca partial waves up 
to l m ax — 10 and p- 48 Ca partial waves up to l ma x = 20 
are included. The maximal total angular momenta con- 
sidered is J m ax = 60. Depending on reaction and energy, 
some of these cutoffs can be reduced without affecting the 
quality of the results. The integrals are discretized using 
Gaussian quadrature rules with about 50 grid points for 
each Jacobi momentum. The Coulomb screening radius 
is Rc — 16 fm. The results are well converged except for 
the breakup cross sections at small angles. In that case 
a special treatment of the pure Coulomb breakup term 
as described in Ref. [l2| improves the convergence. How- 
ever, even with this special treatment including Rc = 50 
fm and proton-nucleus partial waves up to l max = 40 we 
were unable to get better than 20% accuracy for small an- 
gles. We therefore present breakup results neglecting the 
Coulomb force, since the main goal of the present work is 
the comparison of the methods and not the analysis of ex- 
perimental data. The Faddeev-AGS calculations without 
Coulomb force are well converged; in this case l max = 10 
for p-A partial waves is sufficient. 



Label C/ p a U n A nA-bound 



FAGS E d /2 E d /2 


no 


FAGS1 E d /2 E d /2 


yes 


FAGS2 E p E d /2 


yes 



TABLE III: Types of Faddeev-AGS calculations being per- 
formed, the labels used, the energies at which the associated 
interactions were determined and whether a neutron-nucleus 
potential supports a bound state. 

take the nucleon optical potentials U p a and U„a at half 
the deuteron energy in all partial waves and therefore 
produce no bound states in the nucleon-nucleus subsys- 
tems. 

For there to be a transfer channel in Faddeev-AGS, 
the relevant interaction needs to hold a bound state. The 
only difference between FAGS and FAGS1 is the inclusion 
of a neutron-nucleus bound state (by taking a real inter- 
action for the nA subsystem in the partial wave for which 
a bound state exists). Comparing FAGS and FAGS1 for 
elastic and breakup reaction tells us about the impor- 
tance of coupling to the transfer channel. 

In evaluating transfer matrix elements, the calcula- 
tions labeled CDCC take the nucleon-nucleus interac- 
tions (UpA,U n A) in the initial state to be at half the 
deuteron energy Ed/2 and the auxiliary interaction U p b 
to be at the proton energy in the exit channel E p , closely 
mimicking the physical process. Therefore it makes sense 
to consider also the Faddeev-AGS counterpart FAGS2, 
where U p a is determined at E p . We note that FADD2 
and FADD of Refs. IHG3 correspond to our FAGS1 and 
FAGS2, respectively. For completeness we also explore 
the effect of this subtle difference in the initial Hamilto- 
nian in CDCC, by performing CDCC2 (where all proton 
optical potentials arc calculated at E p , including the one 
in the incoming channel). 

A. Elastic Scattering 



IV. RESULTS AND DISCUSSION 

In this section we present results for elastic scattering, 
transfer and breakup. For many of the chosen cases, 
data is available. However, because we neglect spin, the 
direct comparison with data is not appropriate. It is the 
comparison between the methods that is meaningful. 

While the reference calculations are performed with 
the neutron and proton optical potentials calculated at 
half the deuteron energy, there are small useful variations 
in the interactions that should be considered carefully. 
For clarity we describe the various labels that will appear 
in the following subsections in Table Hill For elastic and 
breakup, the standard CDCC and Faddeev-AGS calcula- 
tions are labeled CDCC and FAGS, respectively. These 



The simplest reaction to study is elastic scattering. 
The elastic differential cross section is sensitive to the 
asymptotic behavior of the three-body wavefunction in a 
very particular region of phase space, namely when the 
neutron and proton are bound. It is thus one basic test 
for a complete reaction theory. 

The elastic cross sections for deuterons on 10 Bc, 12 C 
and 48 Ca are shown in Figs. [3J [3] and HI respectively. 
Plotted are the standard CDCC (dashed) and FAGS 
(solid) angular distributions, as well as the results of 
FAGS1 (circles-dotted), when the transfer channel is in- 
cluded. The differential cross section is divided by the 
corresponding Rutherford cross section to allow for in- 
spection at larger angles. Overall CDCC and FAGS 
agree, with exception of d+ 12 C at 12 MeV, where some 
discrepancy is found beyond 40°. The inclusion of the 
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transfer couplings has a minor effect, as can be seen by 
comparing FAGS and FAGS1, which produces only small 
modifications in the distributions, mainly at large angles. 

B. Transfer Reaction 

Contrary to elastic scattering, the transfer process is 
typically sensitive to the three-body scattering wavefunc- 
tion at the surface of the target, and therefore poses dif- 
ferent challenges for reaction theory. Here, all cases refer 
to transfer to the ground state (g.s.) of the final nucleus. 

In Fig. [5j we present the angular distribution for 10 Be 
(d, p) n Be reaction calculated at E d = 21.4 MeV, 40.9 
MeV and 71 MeV, respectively. We first compare the 
standard CDCC (dashed) and standard FAGS1 (solid). 
We observe that at E d = 21.4 MeV, the two methods are 
in perfect agreement, but significant discrepancies ap- 
pear at higher beam energy. To determine how much of 
those discrepancies can be attributed to the choice of the 
energy at which the optical potentials are evaluated, it 
is useful to consider the auxiliary calculations CDCC2 
(dotted diamonds) and FAGS2 (dotted circles). Com- 
paring FAGS2 and CDCC2 (both with U pA at E p in the 
entrance channel) one concludes that the discrepanc y is 
almost independent of the ambiguity in U v a- As in |16f 
we do find an increasing dependence on the choice of the 
optical potentials with beam energy. CDCC calculations 
should be independent of the choice of the auxiliary in- 
teraction. This is true for the lowest energy, but a few 
percent discrepancy starts to appear as the beam energy 
increases. 

In Fig. |6l we show da/dCl for 12 C(d, p) 13 C reaction 
at Ed = 12 and 56 MeV. Just as in the case of 10 Be, 
the CDCC predictions for (d,p) on 12 C at low energy 
provides a very good approximation to the Faddcev solu- 
tion (FAGS1). However, the disagreement at 56 MeV be- 
comes significant (around 20 %). Whereas in 10 Be there 
was a strong dependence of the transfer cross section on 
the choice of the energy at which the optical potentials 
are evaluated, for 12 C, no such dependence exists (com- 
pare FAGS1 and FAGS2 or CDCC and CDCC2) and 
therefore the disagreement is quantitatively robust. 

In FigJ7]we present the angular distributions following 
(d,p) transfer to the ground state of 49 Ca at 56 MeV. 
Small but not negligible discrepancies are found in the 
shapes of da/dCl between FAGS1 and CDCC, accompa- 
nied by a significant dependence on the optical potential 
(FAGS2 and CDCC2), which makes the comparison am- 
biguous. 



We connect the present work with the comparative 
study [HI between the finite-range Adiabatic wave ap- 
proximation (ADWA) method and Faddeev-AGS. For 
that purpose, we include in Figs. |6] and [3 the finite- 
range ADWA predictions (stars-short-dashed). These 
correspond to exactly the same Hamiltonian as the 




20 40 60 80 100120 

(degrees) 



FIG. 2: Elastic cross section for d+ 10 Be: (a) E d = 21.4 MeV, 
(h)E d = 40.9 MeV, and (c) E d = 71 MeV. 
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FIG. 3: Elastic cross section for d+ 12 C: (a) E d = 12 MeV 
and (b) E d = 56 MeV. 
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FIG. 4: Elastic cross section for d+ 48 Ca at E d = 56 MeV. 
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FIG. 5: Angular distribution for 10 Be (d, p) "Be: (a) E d = 
21.4 MeV, (b) E d = 40.9 MeV, and (c) E d = 71 MeV. 
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FIG. 6: Angular distribution for 12 C (d, p) 13 C: (a) E d = 12 
MeV and (b) E d = 56 MeV. 
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FIG. 7: Angular distribution for 48 Ca (d, p) 49 Ca at E d = 56 
MeV. 



target 


E d (MeV) 


Aa-c (%) 


A F _c (%) 


A P _A (%) 




21.4 


-3 


3 ± 1 


6 ± 1 


10 Be 


40.9 


-21 


-36 ± 19 


-19 ± 19 




71.0 


-9 


-53 ± 47 


-48 ± 47 


12 C 


12.0 


8 


6±3 


-2 ±3 




56.0 


13 


-21 ± 1 


-30 ±1 


48 Ca 


56.0 


-5 


39 ± 14 


47 ± 14 



TABLE IV: Percentage differences between the (d,p) cross 
sections at the peak of the angular distribution predicted by 
the various methods, for the three targets here considered, as 
a function of beam energy Ed- 

CDCC calculations performed here, and therefore differ 
from those in [16]. Encouragingly, for the reactions on 
10 Be, ADWA performs just as well or even better than 
CDCC. Fig. [5] is further evidence that, for I = transfer, 
ADWA provides a very good approximation to Faddeev- 
AGS for reactions around 10 McV/u, as concluded in [l6| . 
For 12 C and 48 Ca, ADWA cross sections differ from the 
CDCC predictions by a few percent, usually increasing 
the discrepancy with the FAGS1 results (as compared to 
CDCC). 

Table IIVI provides a quantitative estimate for the dis- 
crepancy between the various methods here considered. 
We provide percentage differences of the differential cross 
sections at the peak of the angular distributions: i) Aa-c 
compares ADWA with the standard CDCC, relative to 
CDCC, ii) A F -c compares FAGS1 with CDCC, relative 
to CDCC, and iii) Ap-A compares FAGS1 with Adia- 
batic, relative to ADWA. The error in Ap-c an d Ap_A 
is estimated from the percentage difference obtained with 
FAGS2 versus FAGS1. 



C. Deuteron Breakup 

In this section, we compare CDCC angular and energy 
distributions following the breakup of the deuteron with 
the Faddeev-AGS counterparts. We consider FAGS and 
FAGS1, i.e., without and with a bound state in the nA 
subsystem. 

Breakup calculations are computationally more de- 
manding than transfer or clastic. As mentioned in 
Sec. IIII B[ the Faddeev-AGS breakup results obtained 
with the present technical implementation are not suf- 
ficiently accurate at forward angles when the Coulomb 
force is included. For this reason, the comparison of 
breakup cross sections is performed switching off the 
Coulomb interaction in both CDCC and Faddeev-AGS. 

In its present implementation, CDCC breakup ob- 
servables are reconstructed from the asymptotic form 
of the three-body wavefunction. This implies that the 
asymptotics of the three-body wavefunction needs to be 
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FIG. 8: Breakup distributions for 10 Be (d, pn) 10 Be reaction 
at: (a) E d = 21 MeV, (b) E d = 40.9 MeV and (c) E d = 
71 MeV. Results for CDCC (hashed bar), FAGS(solid) and 
FAGS1 (circles). 

good not only when the neutron and proton are close 
to each other, but also when they are far apart. This 
poses considerable challenges to our numerical methods. 
For all the CDCC calculations, we went to the numeri- 
cal limit of the methods, mostly dictated by the num- 
ber of partial waves included in the deuteron contin- 
uum. For the nuclear-only CDCC breakup calculations, 
Qmax = Imax = 6 was our model space limit. The er- 
ror bands in the CDCC cross sections are due to the 
truncation of the model space, and have been extrapo- 
lated from the differences found in the cross sections for 
Qmax — Imax — 6 and Qmax — Imax 4. Typically 
the calculations at higher energies have a smaller error 
bar because for those energies the effect of these higher 
multipole couplings is smaller. 

In Fig[5]we present the results for the angular distri- 
bution as a function of the cm. angle of the pn system 
following the breakup on 10 Be at the three energies of 
choice. The hashed band, the solid line and the circles 
correspond to CDCC, FAGS and FAGS1 predictions, re- 
spectively. At the lowest energy, we find that CDCC does 
not reproduce FAGS, even taking into account the error 
estimated by model space truncation. At the higher en- 
ergies, this discrepancy is removed. The insets of FigJS] 
contain the corresponding energy distributions as func- 
tions of the proton- neutron relative energy E pn . Again, 
a very large discrepancy is found at 21.4 MeV while fair 
agreement between CDCC and FAGS is obtained at the 
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FIG. 9: Breakup distributions for 12 C (d, pn) 12 C reaction at: 
(a) E d = 12 MeV and (b) E d = 56 MeV. Results for CDCC 
(hashed bar), FAGS (solid) and FAGS1 (circles). 
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FIG. 10: Breakup distributions for 48 Ca (d, pn) 48 Ca re- 
action at E d = 56 MeV. Results for CDCC (hashed bar), 
FAGS(solid) and FAGS1 (circles). 

higher energies. 

Similar conclusions can be drawn from the comparison 
of breakup angular and energy distributions for reactions 
on 12 C (FigE]). Despite the large error bar in the CDCC 
predictions, there is a striking mismatch between CDCC 
and FAGS in both magnitude and shape of the breakup 
cross sections at 12 MeV. These discrepancies disappear 
at the higher energy. Agreement is obtained between 
CDCC and FAGS for the breakup of deuterons on 48 Ca 
at 56 MeV, as shown in FigfTUl 

The effects of including the nA bound state in the 
transfer channel is shown with FAGS1 (dotted circles). 
By comparing FAGS and FAGS1 we conclude that the 
effects of transfer are not negligible on breakup, particu- 
larly at low energies. 

Although for transfer we found excellent agreement for 
all reactions at low energy (« 10 MeV/u), for breakup it 
is exactly at the low energy that CDCC appears to break 
down. To better understand this, we have investigated 
the relevant contributions to the breakup cross sections 
within FAGS. In particular, we wanted to understand 
the relative importance of the three Faddcev components, 
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namely the deuteron channel, the proton channel and the 
neutron channel, in evaluating the breakup observables. 
We suspected this was the key to understanding the dis- 
crepancies between CDCC and FAGS, because CDCC 
is based on an expansion in the deuteron channel alone. 
Even though in the FAGS calculations there are no bound 
states in the proton and the neutron channels (U n A ancl 
Up a contain an imaginary term), the breakup is still dis- 
tributed between all three Faddccv components. We find 
that, for breakup at low beam energies, the proton and 
neutron Faddeev components are equally important as 
the deuteron component, and there are destructive inter- 
ferences between the various components, at most values 
of Epn, that reduce the cross section significantly. 

At the limiting case where all the energy is in the 
proton-neutron relative motion, the phase space factor 
goes as ~ qA ~ ^/£"p° ax — E pn , where qA is the tar- 
get recoil momentum in the cm. system. This is not 
well reproduced when expanding in terms of a truncated 
deuteron basis, as in CDCC. Given the coordinate choice, 
it is expected that the asymptotic behavior of the CDCC 
wavefunction is at its best when the proton and neutron 
are close to each other, but not in the opposite case. The 
results presented in Figs. [5J HO [TU] arc a manifestation of 
this fact. CDCC breaks down at the lower energy regime, 
where we find broader angular and energy distributions 
and there are larger contributions to the breakup cross 
section from E pn « 

Our previous work [18[ suggested that the discrepan- 
cies between CDCC and Faddeev in breakup arrived as 
a consequence of the existence of a bound state in the 
rearrangement channel. Results shown in FigsJHJ and EO 
demonstrate this is not always the case. Even if there is 
no nA bound state, large discrepancies can occur at low 
energies. The inclusion of a bound state enhances the 
disagreement, as shown by FAGS1 results in Figs. 18191 
and US 



V. SUMMARY AND CONCLUSIONS 

In this study we perform a systematic comparison be- 
tween CDCC and Faddccv- AGS for deuteron induced re- 
actions. The CDCC formalism under study here corre- 
sponds to the standard method introduced in [f|, and is 
based on a truncated expansion on the deuteron contin- 
uum. The Faddccv- AGS method considered in this work 
corresponds to the solution of the AGS equations Q, 
including Coulomb screening and renormalization as in 
[Tlj . For all practical purposes, Faddeev- AGS solutions 
are considered exact and therefore this study should serve 
as a test on the reliability of CDCC. 

We focus on deuteron reactions on 10 Be, 12 C and 48 Ca, 
including a wide range of beam energies. We compute 
elastic scattering, transfer cross sections to the ground 
state of the final system, as well as breakup observables. 
In CDCC, elastic scattering and breakup cross sections 
are obtained directly from the S-matrix, while transfer is 



calculated replacing the exact three-body wavefunction 
by the CDCC wavefunction in the exact post-form T- 
matrix. 

Our CDCC/FAGS comparisons show no immediate 
correlation between elastic, transfer or breakup. In other 
words, finding agreement for the elastic for a given target 
and beam energy does not imply agreement in breakup 
or transfer. Indeed, these processes are sensitive to dif- 
ferent parts of configuration space and therefore, only by 
looking at elastic, transfer and breakup simultaneously, 
can the CDCC method be thoroughly tested. 

Overall, and regardless of the beam energy, CDCC is 
able to provide a good approximation to FAGS for elas- 
tic scattering. The inclusion of a neutron-nucleus bound 
state in the FAGS1 calculations introduces small mod- 
ifications mostly are backward angles. Only for d+ 12 C 
at 12 MeV we found stronger discrepancies in the elastic 
angular distribution between CDCC and Faddccv-AGS. 

The comparison of CDCC and Faddeev- AGS for trans- 
fer cross sections is consistent with the results presented 
in [TU- We found CDCC to be a very good approxima- 
tion of FAGS1 at reactions around 10 McV/u, but not so 
good for larger beam energies. 

As opposed to transfer, breakup observables predicted 
by CDCC are at its best for the higher beam energies 
explored in this work. To reduce the technical challenges 
of the problem, we ignore the Coulomb interaction in 
the breakup comparison. Also, we use exactly the same 
Hamiltonian (CDCC and FAGS) to remove any ambigu- 
ity. Taking into account the estimated error due to the 
truncation of the model space in the CDCC calculations, 
CDCC predictions for the breakup angular and energy 
distributions are in good agreement with FAGS for all 
but the lowest energies considered. For deuteron breakup 
on 10 Bc at 21.4 MeV and on 12 C at 12 MeV, CDCC fails. 
Strong contributions from the proton and neutron Fad- 
deev components, not explicitly included in the CDCC 
expansion, are present when the proton-neutron relative 
energies are large. At low energy, the energy distribu- 
tion is broad, the breakup to scattering states with large 
proton-neutron relative energy is important and there- 
fore CDCC does not perform well. One possible solution 
to this shortcoming is to use the CDCC wavefunction in 
a T-matrix that probes only short distances between the 
proton and neutron, instead of its asymptotic form. 

The present comparisons pose important constraints 
on the validity of CDCC when applied to deuteron in- 
duced reactions. However, it is critical to realize that the 
number of reactions which can be calculated with the 
present implementation of Faddeev- AGS is more limited. 
Indeed, as the mass of the target increases, Coulomb ef- 
fects become stronger requiring larger screening radius 
and more partial waves; this renders the solution of the 
AGS equations in the partial- wave representation impos- 
sible. In addition, the Pade summation technique used 
to solve the AGS equations converges slower at low en- 
ergies (10 McV/u), where couplings are strong, and it is 
not possible to obtain converged solutions even for mass 
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A = 48. It would be desirable to extend present tech- 
niques used in Faddeev-AGS or to develop new methods 
which can overcome these difficulties. 

Finally, one should keep in mind that this work is based 
on a pure three-body formulation of deuteron induced 
reactions. This assumes that the final bound state pop- 
ulated in the (d,p) reaction is of pure single particle na- 
ture. A great challenge and important advance for the 
Faddeev-AGS formulation is the inclusion of degrees of 
freedom of the target, which is necessary for a more re- 
alistic description of the process. 
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